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Abstract 

A new scenario of baryogenesis is presented, within the split fermions framework. Our 
model employs a first order phase transition of the localizer field. The standard model 
(SM), Kobayashi-Maskawa phase induces a sizable CP asymmetry. The usual suppression 
of CP violation which arises in the SM baryogenesis is absent due to the existence of order 
one Yukawa couplings before the fermions are localized in the extra dimension. Models of 
the above type naturally contain B-L violating operators, allowed by the SM symmetries, 
which induce the baryon asymmetry. Our mechanism demonstrates the following concept: 
the flavor puzzle and the SM failure to create the baryon asymmetry are linked and may 
have a common resolution which does not rely on introduction of new CP violating sources. 



1 Introduction 



It was understood long ago that the SM contains all the ingredients required to produce the 
baryon asymmetry of the universe (BAU) [21 

^ ~ 7 x HT 11 , (1) 

s 

where n-Q and s are the baryon and entropy densities respectively. Nevertheless the SM fails to 
explain the observed BAU quantitatively. The lower bound on the Higgs mass 12] implies that 
the EW phase transition (EWPT) is second order (see e.g. [IJ[S] and refs. therein). Thus no 
departure from thermal equilibrium is obtained. Furthermore, the observed quark flavor param- 
eters are small and hierarchical, posing a puzzle, known as the SM flavor puzzle. A measure of 
their smallness can be given by the Jarlskog determinant 0, J = det[Y u Y^, Y d Y ] d ] = O (l(T 19 ) . 
In the SM baryogenesis case, the amount of CP violation (CPV) produced is suppressed by 
J |H1 03 E3 • Thus even under the assumption of a first order phase transition (PT) , the 
possibility that the SM can account for the observed BAU is precluded. 

It is interesting that the failure of SM baryogenesis is directly related to the SM flavor 
puzzle. Recently it was demonstrated in ^2] that a class of Froggatt-Nielsen ^5] models which 
solve the SM flavor puzzle can also account for the BAU. The main idea of P2J is that within the 
Froggatt-Nielsen framework, Yukawa couplings can vary with temperature. Thus it is possible 
that in the early universe the flavor parameters were anarchical, J ~ 1. This implies that CPV 
arising from the Kobayashi-Maskawa (KM) phase, during the electroweak phase transition, is 
unsuppressed. 

In this work we demonstrate that the five dimensional (5D) split fermions framework can 
realize the above idea in a similar manner. In our scenario, fermions are localized due to their 
couplings to an x§ dependent VEV of a bulk scalar denoted as the localizer [141 1151 ITS] . We 
assume that at the critical temperature, two phases coexist where in one of these phases, the 
localizer's VEV is zero (the symmetric phase) while in the other it acquires a non-trivial VEV 
(the broken phase). As in the standard EW baryogenesis, the phase transition occurs through 
bubble nucleation. Outside the expanding bubble wall (in the unbroken phase) the fermions 
wave functions are flat |17l 1181 119j and the Yukawa couplings are of order one. The KM phase 
may therefore be sufficient to explain the observed baryon asymmetry. 

Flavor models typically have two distinct phase transitions. The first is related to the scale 
at which the scalars, which are SM singlets and which control the hierarchy in the flavor sector, 
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acquire VEVs. In our model, we denote this phase transition as the localizer phase transition 
(LPT). The other is the celebrated EWPT. In that context, there are two interesting related 
possibilities: 

(i) The LPT occurs at or below the temperature of the EWPT. 

(ii) The LPT occurs at temperatures well above the temperature of the EWPT |19| . 

Case (i) is very similar to the one discussed in |12| . Case (ii) differs from the SM baryogenesis 
scenario. We show below that it can also naturally avoid all the failures of the SM baryogenesis 
mechanism and account for the observed BAU. In our framework, we do not require new sources 
of CPV or new sources which violate baryon number. The related sources within the SM are 
shown to be sufficient. Furthermore, genericalfy, we expect that in simple split fermions models, 
option (ii) would be realized as follows: naturally the critical temperature for the LPT, Tlpt, 
will be set by the fundamental scales of the problem 0(l/nR) ^ Tlpt ^ (where R is the 
radius the extra dimension). Phenomenologically, the inverse radius, 1/R is constrained to be 
rather high, say above 100 TeV, in order to suppress the contributions to various flavor changing 
processes [201121] . We then expect that Tlpt would be of similar order and therefore much above 
the EW symmetry breaking scale. This implies that the BAU is produced during the LPT, well 
before the EWPT. Let us briefly describe how Sakharov's conditions are fulfilled in our scenario: 
C and CP violation is of order one in the unbroken phase since all the Yukawa couplings are of 
order unity. Out of equilibrium is obtained assuming a first order PT. This can be accounted 
for by, for example, the presence of higher dimension operators in the localizer's potential |19j . 
Baryon violation occurs since there are non-renormalizable operators which are not suppressed 
any further by the wave-functions overlapping. It is remarkable that our mechanism goes 
through even in the most minimal split fermions models a la Kaplan- Tait We shall see 

that, despite the fact that in these models the localizer's couplings are CP conserving |19U22ll2"3j . 
a sizable baryon asymmetry is obtained. 

In section|2]we introduce the split fermions framework and in particular we discuss in some 
detail the specific model we study to demonstrate our mechanism. In section |3] we describe 
our mechanism and explain how we derive a semi-quantitative estimation for the resultant 
asymmetry. In section 0] we show how baryon number violation occurs in our model. This is 
due to higher dimensional B— L violating operators allowed by the SM symmetries. In section [5] 
we present two complimentary methods to estimate the CP asymmetry due to the interaction 
with the bubble wall. Many of the detailed calculations relevant to this section are found in 
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the appendices. We conclude in section H3 



2 The Framework 

In this part we briefly review some of the relevant features of the split fermions framework. To 
demonstrate our mechanism we consider a minimal model with a single localizer scalar. To first 
order the localizer's profile is roughly flat |14[ I15j so that at zero temperature such a model is 
similar to the one in ref j2J with constant odd masses as described below. 

We consider one extra dimension compactified on S 1 /^ with a fundamental domain 
[0, ttR]. The matter content is that of the SM with an additional real scalar "localizer" field <p. 
The relevant part of the 5-dimensional Lagrangian density is given by 

Vj + ^dM^U-V^), (2) 

where M = 0, .., 3, 5, denote SM four-Dirac fermions (which more explicitly are Qi, Ui, di,Li, ei 
with corresponding couplings f%, /y etc., where here and below the representation index is sup- 
pressed), i, j = 1, .., 3 are the flavor indices, M* ~ 10/irR is the fundamental scale in the theory 
and we take 7 5 = — 7°7 1 7 2 7 3 . V((p) is the localizer's potential and a (fi e term is required to ob- 
tain a first order PT ^S] (see also [5] for a recent study of the 4D case) . Due to the 5D Lorentz 
symmetry the /s are hermitian and therefore can all be brought to a basis in which they are real 
and diagonal. Thus CP is not broken by the above action. In more complicated (and realistic) 
models |16] twisting of the fermion wave functions might occur in flavor space [23]) providing 
new CPV sources [191 122j . We shall see below, that even in the limit where these additional 
sources are switched off, the above framework can still account for the observed BAU. 

To avoid the constraints from FCNC |2U1 121j 1/R > 100 TeV is required. Note, however, 
that no upper bound on either 1/R or M* exist. For naturalness we assume that all the 
dimensionless couplings in the theory are of order one. We emphasize that all our computations 
below are done at tree level and we assume that quantum effects are sub-dominant. 

The orbifold boundary conditions are 

<f>(x^, -x 5 ) = -cp{x^,x 5 ), 0(a^, x 5 + ttR) = -<f>(af, x 5 - ttR), 
y(x»,-x 5 ) = ±75^(2^, x 5 ), *(x^,x 5 +ttR) = ±75^(2^, a; 5 -ttR), (3) 

where the +(— ) signs are for SU (2) doublets(singlets). These boundary conditions are intended 
for two things. For the fermions, these conditions project one of the two chiralities of the zero 
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mode, rendering a chiral low-energy theory. For the localizer the conditions impose a non-trivial 
x 5 -dependent VEV, which to our purpose is assumed to be a step function. Thus it dynamically 
plays the role of the odd mass term, a.i of ref j^J (for more precise analysis see e.g. |14M15ll21| ). 
As a consequence of this non-trivial VEV for (j>, the fermions zero-modes are localized in the 
extra dimension, 



where Oj = \fau\/Mj , and fa is the eigenvalue of the Yukawa matrix which corresponds 
to the eigenvector (note that i should not be summed in the above). The small overlap 
between the different zero modes generates the flavor hierarchies. In the model we consider, the 
hierarchies are realized by taking for the weak doublet fields fa > and for the weak singlets 

fa < o EU. 

3 Cosmological Setup 

Below we describe the baryon-asymmetry production mechanism and explain how the above 
framework overcomes the SM difficulties. Our main point here is that, even minimal, split 
fermions models can efficiently produce baryon number through the LPT. It is the fact that 
during the LPT, near the expanding localizer bubble wall, two different phases coexist which 
essentially allows for all of Sakharov's conditions to be realized and lead to a sizable baryon 
production rate. In the unbroken phase CPV and B— L production rates are unsuppressed while 
inside the bubble the relevant processes freeze out due to fermion localization. 

It is important to identify the time, in the history of the Universe, in which the baryon 
production occurs. This is directly linked with Tlpt, the critical temperature of the LPT. As 
mentioned above we naturally expect 0(l/7ri?) <| Tlpt S5 M*. This is also supported by 
the phenomenological requirement of the model as follows. The flavor hierarchy is obtained 
due to the geometrical setup of the fermions in the extra dimension. To account for that, 
the localizer's mass and other fundamental parameters in its potential, should be of order 
O(10/ttR) (to have strong enough localization |15l I16j ). In this case we then expect that the 
critical temperature Tlpt would be of similar order. To avoid constraints from processes which 
induce flavor changing neutral currents the inverse radius is constrained to be rather high, 
1/R > 100 TeV (201 HI]- and therefore the LPT occur well before the EWPT, 




/n>0 
fu<0 



(4) 



Tlpt > ?ewpt ■ 



(5) 
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We shall assume below that ?lpt ~ 1/tvR. This implies that during the LPT, modes with very 
high 5D momentum are Boltzmann suppressed. Hence the dynamics of the baryon production 
can be roughly estimated using our knowledge of 4D field theory. 

Let us discuss in some more detail how Sakharov's conditions are satisfied within our 
framework. 

• CP violation - Outside the bubble wall the quarks wave functions are flat and thus 
the corresponding Yukawa matrices are naturally sizable and anarchical. The presence of 
these sizable, flavor breaking, couplings modify the quark thermal masses |251 EJ 
and imply that CPV is unsuppressed. Since however the LPT occurs before the EWPT it 
is the fact that the quarks have non-universal, CP conserving, couplings to the localizer 
which drives the CP asymmetries. This is discussed in more details below. 

• B violation - Above the EWPT any B+L asymmetry created is washed out by sphalerons. 
It is therefore necessary for B— L violating interactions to exist and contribute to the 
asymmetry. It is remarkable that no new ingredients are required to induce B— L number. 
In fact it was noted a while ago, by Weinberg [25], that non-renormalizable B— L violating 
terms exist in the SM. These interactions are only mildly suppressed outside the bubble 
and are frozen in the bubble after fermions localization has occurred. In fact, the B— L 
operators in our case play the same role which is played by the sphaleron processes in the 
SM baryogenesis. Consequently, we expect that they would yield baryon production in 
our case. 

• Deviation from thermal equilibrium - First order phase transition can be accounted for 
by the presence of a </> 6 term in the 5D effective localizer's potential |19| . This case was 
recently studied in great details for the SM 4D case [3] where it was indeed found that 
such terms induce a first order phase transition and consequently deviation from thermal 
equilibrium. 

Let us now briefly describe how the mechanism of baryon production goes through in our 
model. At Tlpt a phase transition occurs and a bubble of the true vacuum is expanding through 
space. Just like in the SM baryogenesis, as the wall sweeps through space, particles and anti- 
particles hit the bubble wall, interacting with the localizer field. From the four dimensional 
point of view, there is a KK tower of particles and as we show in section the interactions 
with the wall violate KK number conservation, so an even(odd) incoming KK mode is reflected 
into an odd (even) KK mode. CP violation implies that the reflection amplitudes for particles 
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will be different than the ones for anti-particles resulting in an asymmetry in the corresponding 
densities near the interface. This asymmetry induces density of a conserved global charge 
which decays far from the bubble wall due to decoherence effects |27j . The induced charge 
density, in the unbroken phase, biases the B— L violating interactions towards production of 
B— L net density |27[ 17] . This B— L production is shut-off when the expanding bubble wall 
quickly overtakes the region of non-zero charge and the net B— L charge eventually becomes 
our observed BAU. 

An estimate for ub/s may be derived in a similar manner to the one presented in [71 111). 
For simplicity, we assume a 2+1 dimensional problem with (x , x 3 , x 5 ) spanning space. Thus the 
computation of the reflection and transmission of particles in the non-compact space becomes a 
one dimensional problem. Furthermore, we work in the thin wall approximation with low wall 
velocity, v w ~ 0.1. Given the above assumptions one obtains jllj . 



rig / 15 \ f dE 

~ -CB-L 



« + -rC_)-« + -r4_) XA^(£), (6) 



s \ 2ir 2 q* J ^ I 2ttT . 

v a ' n,m 

where ee-L is a suppression arising from the inefficiency of the B-L violating interactions near 
the expanding bubble wall and g* ~ 10 2 is the number of relativistic degrees of freedom. 
is the Fermi-Dirac distribution boosted to the wall's rest frame, n n ± = e^ E ~ VwP ^ + 1 
for a n-KK mode of helicity db in the symmetric (broken) phase. Perturbatively, for a small 
Yukawa coupling to the localizer one can estimate, 



(<+ - n u m _) - « + - n b m _)\ ~ 2n (E) [1 - n (E)] [p^ - p^ j ^ , (7) 

with p^ the momentum in the x 3 direction of the n'th KK mode and uq{E) = {e E / T + 1)~ 1 . 
This is not true in the large coupling limit due to large corrections to the fermion masses, and 
hence to the momentum. Finally, 

Aqp(-E) = Tr I$ nm f{} nm - T{} nm R nm (8) 

is the CP asymmetry related to the reflection coefficients. The quantity R nm is a matrix in 
flavor space that contains the reflection coefficients. For example, i?nm is the coefficient for a 
reflection of an nth KK quark state (with, say, helicity +) and flavor i into an mth KK quark 
state (with helicity — ) and final flavor /. R nm corresponds to the CP-conjugate processes. In 
deriving the equation above, we have used the unitarity relation J2 m (Rn m Rnm + T^ m T nm ) =1 
where T nm is the transmission coefficient for the nth KK mode. As we show in Appendix 
El this relation obtains corrections of order pg JT which are of order unity. Thus we expect 
order one corrections to our estimation of ns/s. In section 15.31 and appendices El and El we 
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present a non-perturbative analysis which supports our claim that the above corrections are 
under control. 

A remark is now in order. The CP asymmetry has contributions from scattering of either 
the quark doublets or the quark singlets. Each of the three fields (Q, u, d) has an independent 
coupling to the localizer, and different Yukawa interactions which produce slightly different 
asymmetries. For simplicity, we calculate below the asymmetry for the quark doublets, however, 
Yukawa interactions mix between the doublets and singlets and since these interactions are in 
equilibrium during baryogenesis, the total CP asymmetry is roughly given by the (weighted) 
average between the CP asymmetries of both doublets and singlets. 

4 Baryon Number Violation 

As we explained in section |3J a successful baryogenesis scenario requires a source of B— L 
violation. Indeed, since the critical temperature Tlpt is well above the EW symmetry breaking 
scale, sphaleron processes are fast both inside and outside the bubble. Consequently, the 
only way to produce baryon number which will not be washed out by sphaleron processes is 
via B— L violating interactions. Furthermore the B— L violating processes which are fast in the 
symmetric phase must be very slow in the broken phase in order to produce the asymmetry. The 
split fermions framework incorporates such B— L violating processes through non-renormalizable 
operators. Below we systematically go through the list of these operators, in order to determine 
which of these is important for our mechanism. 

Starting from the lowest dimensional operator, dimension seven in 5D, one finds the cel- 
ebrated L 2 H 2 . This operator induces Majorana neutrino masses, but is not directly relevant 
to baryon production. Such an operator, however, may not only generate too heavy neutrino 
masses but also washout excess of B— L number and therefore must be suppressed (at least in 
the broken phase) for our scenario to work. This can be realized using additional extra dimen- 
sion, a discrete symmetry (see e.g j^Hj for a recent discussion) or a continuous one ^§1 giving 
the lepton doublets some charge so that this term is highly suppressed. Another alternative is 
if the Higgs is localized away from the lepton doublet. Whatever the mechanism is, we assume 
that it does not interfere with the baryon production mechanism discussed below. 

It was noticed long ago that the symmetries of the SM allow for B— L violating operators _26 
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of dimension 9 and 9^ in 5D. For completeness, we list these operators 

Dimension 9: (dPD M d)(L-f M Q) , (LD M d)(¥-f M Q) , (dP D M d)(e 7 AJ ' d) (9) 
Dimension 9^: (Q E Q)(Ld)H^ , (d?d)(eQ)H^ , (cFd){Lu)H^ 

(d?u)(Ld)rf , (d E d){Ld)H (10) 

where Q(L) stands for a quark(Lepton) doublet, u, d(e) stand for the up and down quark(lepton) 
singlets, H stands for the Higgs field, D stands for a covariant derivative, c denotes charge 
conjugation and M = 0..3, 5 stands for Lorentz indices. In the above we have suppressed color, 
flavor and SU(2) indices. 

Let us estimate the rate for the dimension 9 B— L violating operators of eq. @- Using 
naive dimensional analysis (NDA) we estimate the coefficient of the above operators as (47r) 2 
(a 4-7T per interaction assuming the theory becomes strongly coupled at M*). Thus the relevant 
part of the Lagrangian is given by 



(47T) 
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C = {QD M dd 7 M L + d D M dd^l) + ... , (11) 

where M* is the effective cutoff introduced in section |2J Again, using NDA we estimate the 
rate for B— L violation outside the bubble, Tb-l to be roughly 

(4vr) 4 / T \ 8 

r B - L ~ - — - T. (12) 

In order for the interaction to be in equilibrium we require it to be faster than the expansion 
rate at the critical temperature, 

(4nf /TlptN 8 T l 2 pt 

Tb - l ~— Tlpt {-m7) - Mp7' (13) 

where Mpj is the reduced Planck mass. Using, Tlpt ~ 1/ttR, this yields the following condition 

2 (vriiM,) 7 #1 < 1 . (14) 



(4vr) 3 v Mpj 

The above condition is easily satisfied over a wide range of M* and R. For example for 
RM* = Ap 1 as long as M* < 10 13 GeV the above process is in thermal equilibrium outside 
the bubble. Interestingly, this yields an upper bound for the theory's cutoff and inverse radius 
in our model. The exponentially suppressed overlaps between the lepton and quark fields (and 
between themselves) will ensure that the above processes decouple inside the broken phase. 

Finally, we need to estimate the suppression of the baryon number production which is 
induced by these processes. This suppression occurs due to the fact that the typical time scale 
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related to the bubble wall is expected to be of the order of t w ~ 4/Tlpt w w [3 HZ]. Thus if 
Tb-l^w *C 1 the B— L processes are not efficient enough in producing baryon number near the 
wall, before the corresponding region is overtaken by the bubble. We, therefore, expect that 
the baryon asymmetry to be suppressed by a factor ep-L given by 

e B -L ~ a B -l r B - L t w ~ (-j^r) ~ — ( RM *r 8 > ( 15 ) 

where Ab-l = 2 is the B— L charge difference induced by a single B— L violating process @. 



5 The CP Asymmetry 



The condition for CP conservation in the SM can be formalized by the existence of a matrix 
Kl such that 

k[y u y^k l = (Y u Yj)\ 

K[Y d YjK L = (Y d Yjy. (16) 

In the five dimensional theory at hand, the couplings with the localizer (see eq. (j2J) add the 
condition (considering for a moment only the quark-doublets) 1 

K[f Q K L = r Q . (17) 

As we now show, CP violation in our cosmological context is found (perturbatively) to be 
proportional to Tr[fQ,Y u Yu + l^Yj] 3 . Generically we expect it to be of order one. 

We calculate the CP asymmetry in two different limits as follows: 

• Perturbative three-flavor limit, u/y/M* <C 1/2ttR - where u is the localizer's (roughly 
constant) VEV. This would allow for an analytical computation of the CP asymmetry 
within a full three generation framework. 

• Non-perturbative single flavor limit, u/\fM* <; \/2ttR - an analytic solution is found for 
the reflection amplitude within a single flavor framework. 

The results of our analysis in both limits shows that the CP asymmetry is expected to be sizable. 

Generically, the parameters of our model lie somewhere between the above limits. The exact 

1 In principle this theory contains ten CPV phases |29|. Here we only consider for simplicity the doublet sector 
which contains three such new phase. As shown below our mechanism is successful even if all the new phases 
are switched off which can be achieved naturally as we argue in section HJ 
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result is then understood to have order one corrections, but is not expected to be dramatically 
changed. 



Before moving to the actual analysis let us explain how each of the above limits captures 
different essential ingredients of the dynamics of the interactions with the bubble wall. Such 
an analysis can, in fact, be useful also in other baryogenesis models beyond our framework 
(see e.g [121 1191 !5U] b In order to explain the low energy values of the Yukawa couplings, it is 
required from the quark odd masses to be few times larger than 1 /R . In our framework this 
implies that the localizer's VEV should be roughly u ~ lO/nR^/M* |16| . On the other hand, 
as discussed before, we assume that the critical temperature is of order Tlpt ~ 1/irR. For the 
above value of u, given that the quarks typical energy is of the order of the critical temperature, 
the reflection of the quarks from the localizer wall cannot be treated perturbatively. This implies 
that the non-perturbative treatment can only serve as a rough estimation for the resultant CP 
asymmetry. On the other hand the non-perturbative limit can be solved only within a single 
flavor framework and therefore describes a CP conserving model. It can only be used to estimate 
the reflection/transmission probability and not the reflection asymmetries. 

5.1 The Effective Dirac Equation 

Considering a 2 + 1 dimensional problem, during the phase transition cj) obtains an x 5 dependent 
VEV and the bubble of true vacuum is expanding in the x 3 direction. Hence the relevant part 
of the Lagrangian takes the form 



To demonstrate how the asymmetry is being produced, we consider only the quark doublets, 
namely / = The computation of the CP asymmetry for the other quarks follow in the 
same manner and the total asymmetry is obtained through the thermalization of the doublets 
and the singlets due to Yukawa interactions. Written in matrix form, it is apparent that the 
spin-up and spin-down parts of the Dirac fermion decouple, 



£ 5 = i^a M d M ^i - fijix^x 5 )*^! 



(18) 



where within the thin wall approximation, 




x 3 < 
x 3 > 



(19) 




E + id 3 

f(x 3 ,x 5 ) + d 5 

E-ids 
-/(x 3 ,x 5 ) +d 5 



f(x 3 ,x 5 )-d 5 \ /* + \ 
E-id 3 J \y 3 _J 



/(x 3 ,x 5 )-5 5 \ /¥+ 
E + id 3 J v*i 



) 



(20) 



(21) 
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where the ± subscripts stands for the 4D chirality and we have suppressed flavor indices. 



In order to compute the CP asymmetry, one must take into account fermionic interactions 
with the plasma. To do so, one introduces the notion of quasiparticles (221 HI El EEH , which are 
fermionic excitations in the plasma. As we review in Appendix [XI the effective Dirac equation 
is found to be 

v+ .^ 12 ' t ( *m;*Wv=) = o, (22) 



-^(x 5 )^(x 3 ) + <9 5 V-id 3 + i/2l 
V - l d 3 + i /2l -T(x^)e(x 3 )-d 5 \ 3 5 _ 

-f(x 5 )6(x 3 ) + d 5 r + id 3 + i/2i J x[x > x > - u ' 

where I is the coherence length of the quasiparticles given by / ~ l/<?s^ |31l lllj . Thus the 
quasiparticles are damped as they propagate in space. We also define, 

(24) 

and 

T{x 5 ) = {3u(x 5 )/2Ml /2 )f, (25) 
V = 3(E-Q). (26) 

Here f2 denotes the thermal masses of the quasiparticles which are given, to leading order, by 

nsnanu, 

a% = w! + ^f. (3 + + £ (ny » + ^ , (27) 

nl = + + (29) 

Note that here Y u and Yd are order one couplings as opposed to the four dimensional case. Here 
we have neglected contributions related to the localizer couplings which, although introduce 
further flavor dependence, do not play a role in our scenario. 

We next decompose the fields into KK modes so that the above equations acquire the form 
similar to the 4D case. This is done in detail in appendix We find, 

( V + ids + i/21 -#(x 3 ).F mn -i4 n)N 
I -6(x 3 )F Tmn -p^ V-id 3 + i/2l 



v-id 3 + i/21 -e{x?)T mn - v f>\ , {n)f ^ 

-9(x 3 )F Tmn - pf ] V + id 3 + i/2l 



X (n) (x 3 ) = 0, (30) 
X H (x 3 ) = 0. (31) 
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Here, = n/R, we have suppressed flavor indices, summation over KK indices is implied and 

[ % 2 >i n a fij m-n£2Z + l 

13 (0 Otherwise 

We remind the reader that for quark doublets (singlets), C^?) is projected out. 



5.2 A Perturbative Approach 



We now turn to solve eq. (|3Uj) perturbatively, assuming J rmn is small. This allows us to find a 
quantitative estimate for the size of the induced asymmetry. Since the critical temperature is of 
the order of inverse the 5D radii, it is enough, when calculating the transmission and reflection 
coefficients, to consider only the first few KK modes. All higher modes are unexcited due to a 
strong exponential Boltzmann suppression. In practice we shall only consider the zero and the 
first KK states (note that for each fermion there is a single zero mode and two first KK modes 
with different helicities). 

As is clear from the Green's functions described below, only \ has an incoming zero-mode. 
Hence considering \ and adding a source for the incoming zero-mode, eq. ()3Uj) becomes, 

r ij+ id 3 +i/2i -p^ V^fVA 

-p< n) Vij - id 3 + i/2l J (x 3 )yl " 

f-i5(x 3 )5 ij 5 0m 6 n0 8(x 3 )F% m \ /^ m) (x 3 ) 
V -B(x 3 )J=l nm ) ^("V) 

(n) 

This equation describes free fermions with mass p\ propagating in the unbroken phase, taken 
to be at x 3 < 0. As fermions hit the wall, they interact with the localizer which changes their 
chirality and partially reflect them back. We treat J^ mn perturbatively, describing quasiparticles 
which reflect back and forth inside the broken phase. Due to the non-diagonal form of J :ran in 
KK space, an incoming zero-mode is reflected back into an odd KK-mode (see eq. (|32|l). For 
the reasons mentioned above we consider only an additional first excited KK mode, and so we 
seek an expression for the reflection coefficients, describing the scattering of zero-mode into a 
first KK mode. In order to be sensitive to CP violation, we must take into account multiple 
scattering so as to obtain interference. 
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The free Green's function for this problem is constructed in Appendix [Bj We find, 



G {n) {x 3 -x 




(n) 



(n) 




9(x 3 ' - x 3 )e- ip t x * 



8(x 3 - x 3 ')e 




(34) 



where pg = y V 2 — {p§ ) 2 + i/2l is flavor dependent through "P and is KK-dependent though 

(n) 

p\ . Using the Green's function one easily finds the reflection coefficient. We do this in 
Appendix [U| For now we give the result and note that as opposed to the usual 4D case where 
in order to see CP violation one must expand to fifth order in perturbation theory, here CP 
violation is apparent already at third order. Indeed, 



where here J- = J r01 g s 2 2ttR = (YluRj g 2 s M.J )/ and the sub-indices refer to chirality space. 

Plugging the above reflection coefficient in the CP-odd expression Aq™ given in (|SJ) and 
expanding the Green's function in the flavor dependent part, one finds (see Appendix 



with 5V Q = (T/32)(y u Yj + Y d Yj) and 8V U)d = (T/8)Y Uid Y^ d . This result, as advocated, has 
only mild numerical suppressions arising from the loop corrections to the thermal masses. 

5.3 A Non-perturbative Approach 

The reflection coefficients may be calculated in a very different way. Instead of the perturbative 
approach which is sensitive to the CPV in the theory we can solve a simpler reflection problem 
semi-analytically. This is important since in order to solve the flavor puzzle we require / of 
order a few which is beyond the perturbative regime. The problem with the following approach 
is that the calculation can only be done in a single flavor framework so CP violation is absent. 
Our attempt here is therefore only to have a quantitative estimate for the magnitude of the 
reflection coefficients. We then assume that when promoting the problem to accommodate three 
generations, the reflection coefficients become complex while their magnitudes only change 




(35) 




-2 



(36) 
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Figure 1: Estimation of A(E) 01 , the probability of a zero mode to be reflected into a positive 
helicity state, in a single flavor model for 1 < ER < 2. 

by order one coefficients. The explicit calculation is done in Appendix [D] and here we only 
summarize our results. 

For the calculation we adopt the following strategy. We apply the thin wall approximation 
and solve the Dirac equation exactly inside and outside the bubble for the whole KK tower. 
Since the dominant contributions come from the lowest KK modes, we consider only the zeroth 
and first KK modes as on shell modes (i.e. only these modes will contribute to the reflection 
and transmission asymmetry). We then match the wave functions of the incoming zero mode 
and the first KK mode outside the bubble to the outgoing ones inside the bubble. This allows 
us to extract the reflection coefficients. While this procedure is clear and analytic, the final 
stage of the actual matching involves the matching of a single state outside the bubble to an 
infinite tower of (virtual) states inside the bubble which can only be done numerically. We 
verified our numerical results in several limiting cases as described in the Appendix. Our main 
result is that a sizable reflection is found over a wide range of energies |TJ||. This is shown in 
fig. ^ which shows the reflection coefficient for an incoming zero mode. It is apparent from the 
figure, that the reflection coefficients vanish as the energy drops below 1/R. This is expected 
since at such energies only the zero-modes are excited, and thus cannot be reflected by the wall. 
In our analysis we only consider the first KK states so that it applies only for w < 2/R . 

As a final remark, we list the correct unitarity relations between the reflection and trans- 
mission coefficients, which are calculated in Appendix lEl 



r fc |i? 01+ | 2 + |7oo| 2 = l, | J Ri-i+| 2 + -|7i- | 2 = l, -|i? 1+0 | 2 + |iW| 2 = l 



(37) 
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where = w 1 — E i R z and we have explicitly written the helicities of the incoming and outgoing 
states. As was remarked in section |3 these exact relation would imply that the final expression 
for ub/s is slightly altered when the realistic values are considered for the model's fundamental 
parameters. Nevertheless our-order-of-magnitude estimate for the asymmetry is insensitive for 
such order one effects. 



6 Discussion Sz. Conclusions 



We begin this part by presenting our final result. In section |1] we showed that, ee-L, the 
suppression from the inefficiency of the B— L production near the bubble wall is given by 
^B-L ~ l/[^w(-R-^*) 8 ]- in section|S]we find that, A^p, the CP asymmetry in the reflection rates 
is mildly suppressed and is given by Aqp ~ 10~ 2 . Boltzmann suppression and other numerical 
factors provide another suppression of order O(10~ 1 ). We therefore estimate the ratio between 
the baryon asymmetry and entropy to be given by 

^ ~ — e B - L ~ 10- 10 x (Ac x RM*)~ 8 . (38) 
s g* 

where g* ~ 10 2 and Ac ~ 0.22 is the Cabibbo mixing angle and for comparison with Froggatt- 
Nielsen type models it is convenient to consider l/i?M* in units of Ac - This is since in the 
Froggatt-Nielsen framework Ac sets the ratio between the Froggatt-Nielsen scalar VEV and 
the effective cut-off scale. Thus our estimation for the baryon asymmetry is consistent with the 
observed value. 

The flavor puzzle of the standard model (SM) is related to the fact that most of the flavor 
parameters are small and hierarchical which hints towards physics beyond the SM. Several kinds 
of solutions have been suggested to solve the above puzzle. They all require additional fields 
with new dynamics. Our main message in this work can be stated as follows: In a cosmological 
context the additional dynamics, which exists in flavor models, naturally incorporates order one 
Yukawa couplings in the early universe. This could result with a baryon asymmetry produced 
with the KM phase as the only source of CPV. Thus it is plausible that the flavor puzzle and 
the SM failure to create the baryon asymmetry are linked and have a common resolution which 
does not require the presence of new CP violating sources. This is not to say that no new CP 
violating sources exist beyond the SM. We merely point out that the KM phase may play an 
important role in the production of the baryon asymmetry. 

In ^2] the above concept was introduced and realized in Froggatt-Nielsen type models in 
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four dimensions. In this paper we illustrate how a similar behaviour is obtained within a different 
framework which solves the flavor puzzle, that is the split fermion framework. Assuming a first 
order phase transition in which a field (the localizer) responsible for the localization obtains a 
VEV, a net baryon number is produced. CPV occurs through interactions of fermions with the 
bubble wall. Furthermore non-renormalizable B— L violating interactions, which are allowed 
by the SM symmetries, induce baryon production outside the bubble wall. These interactions 
decouple in the broken phase due to fermions localization and thus play a similar role to the 
sphaleron processes in the SM case. 

We argued above that no new sources of CP violation are required for our mechanism 
to work. In our model CP violation was found to be proportional to (at the perturbative 
level) Jf = Tr(/)det[/, Y'Y] where / and Y denote the Yukawa couplings of the fermions to 
the localizer and the Higgs respectively. Both / and Y are matrices in flavor space and are 
expected to be anarchical (we omit the flavor and representation indices for simplicity). Let 
us demonstrate that indeed one can construct a natural model in which the KM phase is the 
only CP violating source. Consider a model in which the two Yukawa matrices are promoted to 
be fields which transform as bi-fundamentals of the corresponding SM U (3)q x U(3) u x U(3)d 
flavor group. In the case in which only these two bi- fundamental fields break the flavor group 
the non-universal part of / will be given, to leading order, by a bilinear combination of the 
Yukawas. Thus it is clear that any CP violating quantity in the above theory (in particular Jf) 
must be proportional to the commutator of the two bi-fundamental fields. Thus the only CP 
violating source is the KM phase. This should hold even in the presence of quantum correction 
and therefore technically natural. 

In |19j it was demonstrated that realistic split fermion models may induce lepton production 
from different CP violating sources [22] , namely through twisting of the fermion wave function 
in flavor space [23] • It is interesting that the effect described above is complimentary to the 
one induced via twisting and in general both are expected to be present. 

Our mechanism requires that the critical temperature should be of the order of the inverse 
compactification scale, 1/R. Phenomenological constraints typically yield 1/R > 100 TeV. 
This implies that our framework allows for a rather low reheating temperature, after infla- 
tion, solving the various moduli/gravitino problems. Finally, we note that our model which 
is based on Kaplan- Tait type models |21j . with a constant localizer VEV, is compatible with 
five dimensional supersymmetric split fermion models (see e.g. \3'2\). Hence it can easily be 
supersymmetrized. 
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APPENDIX 

A The Effective Dirac Equation and KK Decomposition 

As explained in jllj . the dispersion relation for quasiparticles may be found by looking for the 
poles of the full propagator, including the self-eneries. Taking thermal self-energy of the form, 

E(E,p)= 1 °a(E,p)-b(E,p)i-p, (39) 

one is then looking for a solution to 

det[Sv 1 -Z(E,p)]=0, 

with Sq 1 = j°E — 7 • p. Considering the 2 + 1 problem at hand, the 

E = a(E,p) ± [1 - b(E,p)] y/pl+pl 

In the four dimensional case, the self-energy can be linearized as 

E{E,p)~j°(2Q-i/3l-E)-7-p/3, (42) 

where as mentioned in section |5~T1 f2 denotes the thermal mass of the quasiparticle and I is the 
coherence length of the quasiparticles given by I = 1.1/ 'g 2 s T [3D EH- This is of 

course also true 

in the broken phase where the theory is effectively 4D. In the symmetric phase, the theory is 
effectively 5D and thus we expect corrections to the above. Nevertheless, at the temperature 
of phase transition, the 5D linearized self-energies are expected to be affected only by a few 



(40) 

solution to the above is 
(41) 
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higher KK states and are thus expected, at zeroth order, to admit the same behaviour. With 
the above, one obtains the dispersion relation, 

% 1 



E = n--±-y/pl+pl (43) 

The effective free Dirac equation is therefore obtained by taking — ► "f^p^ — S which has 
the effect of taking E — > 2{E — £1 + i/6l) and p — > 2p/3. Thus, plugging these in eq. ()2Uj) and 
massaging a little one finds the effective Dirac equation given in (|22 |) . H23|) . 

Next we KK decompose the Dirac fermions. The KK reduction of the wave function is, 

00 



^ ± (* 3 ,* 5 ) = ^uf n (x 5 )^ n J(x 3 ), (44) 



"ml- 1 

n=0 

^5 



where the ± are the two chiralities of the Dirac fermion, and uf n = af n e ipbX + bf n e ip5X . The 
boundary conditions for uf n are given by (again, considering, for example, the quark doublets), 

uf n (x 5 = + ,7rR-) = ±uf n (x 5 = 0-,7TR + ), (45) 
d 5 u+\ x 5 =0+ ^ R - = f^x 5 = + ,7rR~)u+ n (x 5 = 0+,ttR~), (46) 
d->u~ n \ x s =Q +^ R - = d 5 u~ n \ x 5 =0 - i7rR +. (47) 

Hence in the unbroken phase, 

2 \ 1/2 /nx 5 \ ( 2 \ 1/2 /nx 5 \ . , 

^r) cos (^J' u -={^r) sin Uj- (48) 

It is important to note that under this decomposition, the thermal masses are KK diagonal. 
To see this simply note that we assume the Higgs field is x 5 independent. Thus the Yukawa 
interactions decompose into 

Y-jHQd — YijH (O?^ + Q^Sf) . (49) 

n 

Here we replaced (temporarily) by Q and d to stress that the interaction is between two 
distinct fields. Note that Q and d have orbifold boundary conditions with opposite sign. The 

(n) 

above in particular means that f2± are independent of n. Moreover, Q + = fi_. 



On the other hand, the interactions with the localizer are not KK diagonal. We have, 



urn 



with T given in eqs. I|25jl. (|3*2")l . In the above, the transpose is with respect to KK number. 
Note further, that due to the complex conjugate term in the Lagrangian, T can be taken to be 
Hermitian with respect to flavor. 
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B Computation of the Green's Function 



To construct the (free) Green's functions we solve: 

G§V-z 3 ') = iS{x a - x^Stk, (.11.) 



(T + )ij+id 3 -pf> \A(n),ja ,,-./, _ :x , ,;■ 



- id 3 -P5 



(n) 



G^(x 3 -x 3 ') = i5(x 3 - x 3 ')6 lk , (52) 



where G^(x 3 — x 3 ') is a 2 x 2 matrix. The solution is found by a double Fourier transform. 
Defining, 

G$(x 3 - x 3 ') = j dke ik ^ 3 -^G^(k) (53) 
G$(x 3 - x 3 ') = j dke ik ^- x ^G^(k) (54) 
we obtain (using V+ = V- = V), 

6ft (k) = 27ri( r - { * -*f]) = -2ni | -^r| (55) 



/ jk v . „ . i 



= « P| = -2* 





fc 2 


-fcn 2 




-k+r 


fc 2 


M n) ) 2 






fc 2 


-(pi n) ) 2 




k+r 


fc 2 


-(?4 n) ) 2 



jk 



jk 

where (p^) 2 = V 2 — (p^) 2 - This is precisely the Green's function of a massive fermion, 
(7 M p M + m)/(p 2 — m 2 ) with m = p^ . The inverse Fourier transform depends on the contour 
chosen for the integral which is equivalent to choosing a boundary condition. There are four 

in) 

different contours going above or below each of the two poles k = ±p 3 . For reasons to be 
understood below, we choose the contour which goes above k = —p^ and below k = +p^ ■ 
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After a Cauchy integral one obtains, 



x 3> ) = 



(V+p^ 


(n) 

Ps 


7 \ 

2p"> 


/ \ 


Ps > 




\ 2p<"> 




/v-pP 


P w 




2p<"> 




7>+p<" 




7 \ 


(V-pf 




2p<"> 


2p<"> 


pP 


7>+p<" 


\ 2 P ^ 


2p<"> 


/V+p\ n) 


p<"> 


I 2p^ 


2pW 


1 


7>-p(" 


V 2p<"> 


2p(") 



6(x 3 - x 3 ')e ip z 



9(x 3 ' -x 3 )e-^ )x3 , 



jk 



6(x 3 - x 3 ')e ip s 



.(n) x 3 



jk 



8{x 3 ' -x 3 )e- ip i n)x3 . 



jk 



(57) 



(58) 



The above matrices can be diagonalized. One finds one left-moving mode and one right- 
moving mode. The diagonal Green's functions are given by 



GgV-* 3/ ) 



<^V-* 3/ ) 



^9{X 6 - X 3 '^ 







V 

/ V9{x 3 ' - x 3 )e-^ n) - 3 
{ 





v 8(x 3 ' - x 3 )e"^ n) ^ 3 



Pa 



v 0{x 3 _ ^lypf^ 3 

Pa 



jk 



jk 



and are obtained through a change of basis with the diagonalizing matrices, 



» 



U^=9{x 3 -x 3 ') [ P3 +V 

c/W = e(x 3 ' - x 3 



{n) 

p\ 



(n) 

1 

„(") 



+ 9(x 3 ' - x 3 ) 



(n) . ^ p^-P 

Ps 

vf i 



+ 9{x 3 - x 



3/\ 



(n) 
^5 

Ps 

P\ - 1 



(n) 

Ps 

-1 



(59) 



(60) 



(61) 
(62) 



Note that with this choice there is a well defined p§ — > limit, but the matrices are not 
orthogonal. One can choose an orthogonal matrix by normalizing the eigenverctors. Note also, 
that there is only an incoming (outgoing) zero-mode for x(x) since there is no X-\ X-^- Finally, 
it is clear that the Green's functions fulfill, 



G+(-oo) = G_(+oo) = 0. 



(63) 
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C A Perturbative Computation of the CP Asymmetry 



Our goal is a perturbative solution of the reflection coefficients for the scattering described by 
eq. As we argued in section IS"2l the lowest order su(3)-invariant expression is of the 

form Tr(j4 3 ) = 3det(j4) for A E su(3). Thus the first term to contribute to the asymmetry, 
A C p ~ i*oi-Roi, is of form Ti([T,Y^Y} 3 ) and can therefore arise from fourth order in T . The 
reflection coefficients, should therefore be calculated to third order in T . 



Using the Green's functions, one transforms eq. (j33j) into an intergral one, 



(64) 



F nm \ l*i im) (x 3 ^ 



where we suppressed the flavor indices. 
The solution is found iteratively, 

\9£ n \a?)) V J 



(66) 



, T ,l(n)/ 3\\ 2 r- / fl r«m\ 



x 3 ' - x 3 ") 



(68) 



-^ m0 o J { 



G^(x 3 ' - x 3 ") (_^ Tmi ^f) G«(x 3 " - x 3 "') ( J TW 

\ jk / \ jk 



Note that as opposed to the 4D case, here, in general, all orders contribute to the reflection 
coefficients due to the non-diagonal form of the Green's functions. Nevertheless, the special 
form of J rmn J eq. (|36|h permits only odd (even) non- vanishing orders in the expansion for 
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an odd(even) outgoing KK mode. Even with this simplification, the general expressions for 
the reflection coefficients are very awkward. To simplify, we consider only the zeroth and 
first KK modes. As explained in section 15.21 this is justified since these two modes have the 
lowest Boltzman suppression factors. Furthermore, numerically, J :nm is larger for n, m = 0, 1. 
Focusing on these KK modes, we see that the integrals above are highly damped for x 3 jl S> 1 
due to decoherence. We thus estimate them by replacing x 3 with the decoherence length 
/ = g~ 2 2trR. Doing so results with eq. ()35|). 

The only flavor dependence in the Green's functions comes from Y'Y in V (by Y^Y we 
mean Y^Y U + Y}Y d for the SU(2) doublets and Y u{d) Y^ d) for the singlets). Using V ~ V d + 6V 

(") - . / ( T) , _i_ KT>\2 _ J n ) 



and pf } = J(V d + 5V) 2 - pf> + i/2l with bV oc Y+Y we expand 



G (n) =G ^a a (l + J2C%l(l5V) k ) (70) 

k=l 



where, G\j 10,9 is the flavor diagonal Green's function obtained by setting Y^Y to zero, and for 

(n) 

the sake of brevity, we denote the flavor independent coefficients by C\- k . Thus for example, 
C[% = i k /kl etc. 

Below we choose the flavor basis such that T is diagonal. Note that the diagonalizing 
matrix does not contain (by definition) any CP phase, namely, T is real. So to order J 74 , 

R t o _ r m T ( r w\ r w , r w\ r w \ Tr {v) Tr m Tr (p) 

-Koi^Ol — <^ii -r I (-x 12 <j 12 + Lr 22 Cx 22 I 7-U 22 7-Or n 

(0)t xv^Wt -r^(°)t t I V^Mf nW , nW\ -crffl 



' ^11 ^^22 ^ \ U 12 U 12 "+* < - T 22 ^22 J • ALr ll V ^ 

Each term in eq. (|71|) contributes to the CP asymmetry seperately. In expanding the 
Green's functions, the first term to contribute is of order 5V 3 . Using the facts that SV is 
hermitian and that -Roi is obtained from i?oi by replacing 5V with 5V T we find, as an example, 
for the first term in the first line, 



A C p oc 



r (0)* ( r (0) r (i) r (0) r {\)* _ r (i)* r {\) _ r {i) _ r (i)* 

°11,1 I °ll,l L/ 12,l "r" U 11,1 L/ 12,1 u 12,l°12,l L '12,2 u 12,2 



+ c ii*2 ( C u,i + C i2,i)\ Tr [pV^TbVT 3 - 5V 2 F 3 5VF) 

oc Tr 5V] 3 ) = Tr (JF) det ([JF, 5V\) . (72) 

where the proportionality constant is just the flavor diagonal product of Green's functions. 
Similar contributions arise from the other terms in ()71|) and all in all we obtain the expected 
result A C p oc Tr(.F) det([J r , 5V]). 
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D Non-Perturbative Calculations of the Reflection and Trans- 
mission Coefficients 



In this part we calculate the reflection and transmission coefficients of the fermions in the non- 
perturbative regime (large Yukawa couplings). This is possible at the expense of working within 
a single generation framework. Thus CPV is absent and therefore thermal effects, in particular 
thermal masses, are neglected for simplicity. Our rationale here is that our results below, for 
the reflection coefficients in the single flavor case, would also hold when promoting the model 
to three generations ^H]- The main difference is the appearance of CPV phases in the thermal 
masses [see e.g. eq. (|27j)]. These are expected to be only mildly suppressed by an order one 
factor [1/8(1/32) for the quark singlets(doublets) in the 4D SM case] which we shall take into 
account in our final estimation for the resultant BAU 2 . Consequently, below we shall assume 
the following dispersion relations 



where here and below, in order not to confuse with the perturbative approach, we use m n = p 5 . 

As before, for X3 < the system is in the unbroken phase while for X3 > X3 the system is in 
the broken phase. In order to calculate the reflection and transmission coefficients we solve the 
Dirac equations, for the two regions separately and sew the solutions at X3 = 0. As mentioned 
before the spin up and down states do not mix in the presence of the localizer's VEV and thus 
below we only consider the wave function for the spin down states where similar derivation can 
be applied to the spin up case. 

D.l Unbroken Phase 

For x 3 < 0, the relevant positive frequency solutions are described by plane waves propagating 

in the positive x 3 direction with momenta k^. Since we are dealing with T ~ R^ 1 we will mostly 

be interested in the dynamic of the lowest KK states. For each representation this implies a 

single zero mode which can be either incoming or outgoing, but not both. For simplicity we 

only consider particle with energies less than 2/R and therefore in addition to the zero-mode, a 

massive Dirac state corresponding to the first KK level can be on shell (while higher KK modes 

are Boltzmann suppressed). Consequently we also consider incoming or outgoing massive states. 

2 Note that in |19|. CP violation originated from a different source, namely twisting 1221 1231 . so that the above 
mild suppression is absent and an asymmetry was induced even in the absence of these thermal effects. 
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From (|18H26l44"ll48j) we find that the solution is characterized by a KK number n (consid- 
ering as an example an incoming spin-up fermion): 



x (x 3 < 0,n,k 3 ) = APv( fc 3z 3 -i?t) 



cos ( n%r ) yJE + k 3 



sin ( n^r ) \fE~- 



where is a normalization constants and 

A, 



E = J m 2 n + k\ , 



n 

R? 



(74) 



(75) 



Similar expression would hold for the outgoing (reflected) particles which can be obtained from 
the above via changing k 3 to — k 3 . As before, we see that the reflected states do not contain 
a zero mode for the spin-up fermions. Similarly, for spin-down fermions there is no incoming 
zero mode. 

We can now extract the reflection coefficients at T ~ 1/R by considering an incoming zero 
mode with E > 1/R (to allow for an on shell reflected state) which is reflected to a KK tower 
of massive states. Thus in the unbroken phase the wave function of the fermions is given by 



^ (x 3 < 0) 



-iE t 



cos ( 2L- 



8111 



COS 



8111 



ttR 

E-k\ 
ttR 



R 



E-k% 
ttR 

nR 



+ he 



Ex 3 



- 1 

ttR VO 



(76) 



where above and below we use normalized wave functions for the part which depends on a; 5 , we 
set the normalization constant for the incoming particle to one and unless otherwise specified 
the summation of n is implied and taken from one to the highest KK mode. As mentioned 
above we only consider incoming modes with E < 2/R. Thus one can distinguish between 
two cases: (i) the incoming mode is a zero mode, which implies Iq = 1 and I\ = 0, or (ii) the 
incoming mode is a first KK mode which implies setting I\ = 1 and Iq = 0. 

D.2 Broken Phase 



For x 3 > the general solution depends on f(x 3 ,x 5 ) and cannot be found analytically. Recall 
that, for simplicity, we adopt the thin wall approximation (in the x 3 direction) and assume a 
step function (in the x 5 direction) for the localizer, eq. ()19|) . To simplify our notations below 



we use: 



fij {x , X 





/ 



x 3 < 
x 3 > 



(77) 
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Thus, for a given KK mode, the positive frequency solution is given by: 



X 



where in this case /c 3 2 , the square of the momentum in the x 6 direction, can be either negative or 
positive, corresponding to on shell and virtual modes respectively. Solving for the x 5 dependent 
part, using eqs. l|18l2fi|44T48JI . we find that there are two types of solutions: an exponential 
solution corresponding to a zero mode, 



( x 3 > _ N b J{k'^-Et) ( <(^)V E + k T 3 



(78) 



u+(x 5 ) 



Jx 5 



u (x 5 ) = 0. 



and an oscillatory solution which corresponds to higher KK modes 



E 2 — k' 3 — fc| 



(79) 



u„ \ x" 



n- 



irR(f 2 R 2 + n 2 ) 



cos (k£x 5 ) + ^ sin (££x 5 ) 



E 2 = k' 2 + kf 



K 2 + r + ^- 



n 



(80) 



— • sin (kl-x ) , ^ — r 3 -r re 5 — r 3 , ., , R2 

Hence in the broken phase the wave function of an outgoing fermion with positive spin is given 
by 



X {x 3 > 0) 



-iEt 



e 3 



( " K I^J+]4 sin (^ 



sin 



i 2 +f 2 R 2 



E-k' 3 n 
ttR 



(81) 



+ T e 



iEx 3 +fx i 



2Ef 



e 2nfR _ 1 V 



where 



, 3 - v ^--m n - ^ , R2 . 

It is important to note that since in our framework fR is large (to produce significant local- 
ization), fR ~ 5 — 10 , the non-zero modes in the broken phase are very heavy. Thus only the 
zero mode is produced on-shell and be transmitted through from the unbroken phase. 



E 2_ f 2 



n- 



(82) 



D.3 Computation of the Reflection and Transmission Coefficients 

We can equate the expressions for the wave functions in the broken and unbroken phase, and 
extract the reflection and transmission coefficients. Below we distinguish between the two 
possibilities for an incoming state, namely between an incoming zero mode and an incoming 
first KK mode. As a consistency check we also verify below that our solutions make sense in 
the / — > limit which can be computed analytically. 
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D.3.1 Reflection Coefficients for an Incoming Zero Mode 



The reflection and transmission coefficients are found by requiring that the wave function is 
continuous at x s = 0. This is done by comparing the expressions in Q76|) and (|81[). We note 
that the functional dependence on x 5 for the lower component (negative helicity) in both cases 
is equal up to a kinematic factor. This implies the following relation between the coefficients: 



Rl 



E + k% 



E + JE 2 



E 



El ~ P 



(83) 



Substituting the above relation into the upper component of ()76l81j) we get the following equa- 



tion for the reflection coefficients (after multiplying both sides by \j '^): 



1 + 22 R n e n cos 



nw 



Y 

,2nY _ ^ + 



^ ^ Rn&n-A n 



Y 



cos (nw) H sin (nw) 

n 



(84) 



where 



w 



x 5 /R, 



X = ER , Y = fR 



A,. 



lE + k% E + k' 3 n 



Y 2 + n 2 



E-kl 



E-k' 3 n ' 



n^O. (85) 



In order to find the values of the various coefficients we expand all functions in the cosine basis, 

- Ylx 1 2Y ^ (-l) n e Y7r - 1 
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1 e 1 

7T 



Y 
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7T ^ 



Y 2 + n 2 



cos (nw) , 



sin (nw) 



1 l-(-l) n 2n 1 - (-l) n+m , v 

H > ^ ^ — cos (mw) , m f= n 



7T 



)1 



— m^ 



(86) 



where in the second line if m = n the coefficient is zero. Substitute (|86[1 into (|84|) and collecting 
terms we get: 
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ttY 
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e YlT _|_ I ^ L , n 
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2Y 1 - (-\) n + m 
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(87) 
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where summation over n, m begins at one. 



We can use the relation for the zeroth element to express Tq as a function all the other 
unknowns, 



'Yvr e Y7T + 1 
~2~ ' e Y7T - 1 



Y 



1 y R n e n A r , 

IT ^— ' 



l-(-l) 



rr 



Substituting (|85|) back into (|8T)) we get, 

nm ^ 
1 - 
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7rn 



Y 2 +m 2 
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y 2 + m 2 



cos (mw) . 



Let us write the above relation in a matrix form, describing the set of equations for the various 
Fourier coefficients 



M Rn 7? — V 



(90) 



where M^ n and V n can be directly extracted from (|89|) . Inverting M yields the solution for the 
reflection coefficients 



R n = {M±) l V n , 



(91) 



This of course is done numerically. 



D.3.2 Reflection Coefficients for an Incoming First KK Mode with Spin Up 



As before, we write the relation between the reflection and transmission coefficients, this time 
for higher KK modes, 



E +kV 



Eo-k 3 

Again we substitute the above to get, 
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(93) 



30 



where 



B r , 



Y 2 + n 2 



! E + k' 3 n 



(94) 



Expanding and collecting terms, using l|86|). one obtains, 
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which leads to, 



! Yn e Y7r + l 



2 e y7r - 1 



2Y 



7T 



-ei 



y i — f— i) n 

Si + - V i^e^^ • 

7r n z 



Substituting (jHEJ in (jMI) we get, 
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which may again be used to extract the reflection coefficients. 
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D.3.3 Reflection Coefficients for an Incoming First KK Mode with Spin Down 



This time, 




/ " V ^"+ft; n ¥ zl , 

n = 1 ; 



31 



which, after substitution, gives, 
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Finally, repeating the same procedure as before, one gets 
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D.3.4 Check Y -> 



In order to obtain some intuition and check our results, let us check the above in the small 
energy, small barrier limit (see (|85|0. 

Y -» 0. (102) 
In that limit, to leading order, relation Q89|) is simplified 



2Y i — (_i) m 

i2 n e n (^4 n - 1) 5 mn cos (rmo) = — 5 cos (mw) , (103) 



and we find the following expressions for A n (|85|). 



A n = -l + 2-J^-l + 2^ + 0(y 2 ) . (104) 
Under the above approximation the explicit matrix M mn ()90|) can easily be inverted and 
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we find the following approximated solution for the reflection coefficients: 3 

AY 

Rln-l = — -j , R2n = 0. (105) 

vr (Mn-i - 1) e 2n -i (2n - 1) 
Furthermore, using (|88|) we also have the transmission coefficient, 



T Q = [l + 0(Y 2 )] 



8Y 2 ^ A 



2n-l 

2n 



' x (A 2 „_i - 1) (2n - l) 4 



(106) 



The above simplifications and the generalizations to the reflections of other modes, allow 
us to apply several checks for our numerical results. 

• Vanishing of odd reflection coefficients for small Y. Verified numerically. 

• When X < 1 we expect iV^ to be one. Verified numerically. 

• For 2 > X > 1 and 7 <C 1 we expect N^ to deviate from unity only at the order of Y 2 . 
The reflection coefficient of the first KK mode, is expected to deviate at order Y. Verified 
numerically. 

Another way to verify that our numerical results is correct is to check that it is satisfying the 
unitarity condition as we discuss next. 



E Unitarity Conditions 

We now wish to establish the unitarity conditions for our five dimensional scattering process. 
To obtain intuition, we consider the relation between the reflection and transmission coefficient 
in the following two cases: (i) a non-relativistic scalar - one particle quantum mechanics de- 
scribed via Schroedinger equation and (ii) relativistic spin half fermion - one particle quantum 
mechanics described via Dirac equation. 

E.l Schroedinger Equation 

We wish to calculate the reflection of transmission coefficients of a particle in a step function 
potential: 

V(x) = V6(x 3 ) . (107) 

The reason that the even coefficients vanish is probably related to fact that the potentiai is symmetric with 
respect to the reflection around the point w = n/2. Thus to first order oniy transitions between states that 
respect this symmetry are aiiowed. 
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The Schroedinger equation reads: 

(-V 2 + V(x)) ip = Eip. 
For an incoming particle with momentum k = (0, 0, k) the solution is given by 



^(x 3 ) = 6(-x 3 ) 



e ikx3 + re- ikx3 



3\+Jk'x s 



+ 6(x i )te 



k' 2 = E- V. 



Continuity of the wave function and its derivative requires, 
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To obtain the unitarity condition, we consider the Schroedinger conserved current equation 
(for static density all over space) which is given by 



(,* \ |0O ri 

-co — u ' 



v • (v>*vv> - ^vf) = o. 

Integrating the above from — oo to oo and omitting the interference term we find 

V*w - vv^* 

which implies for (|109j) the relation: 

\r\ 2 + R k \t\ 2 = 1, 
and clearly ()110|) is consistent with the above. 
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E.2 Dirac Equation - 4D Case 



Let us first consider a reflection problem from a one dimensional wall embedded in four dimen- 
sions. The Dirac equation reads, 



7° [V • 7 + V(x)} i\) = Ei;. 



(114) 



For a massless incoming spin down particle with momentum = (k, 0, 0, k) the solution is (in 
the chiral basis), 
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Note that in this case (which is different from the 5D case) the transmitted particle is massive 
while the reflected one is massless. From continuouity of the wave function one finds: 



r = 



l-Rk 
l + Rk 



l + Rk 



k' 



(116) 



Again, unitarity condition is found through the Dirac conserved current equation, 

V • (^) = 0, (117) 

which gives 

(«) l-oo = 0. (118) 



We thus obtain the condition, 

\r\ 2 + R k \t\ 2 = 1. 
E.3 Dirac Equation - 5D Case 



(119) 



We now move on to discuss the 5D problem. For finite temperature of order (1/R) we assume 
it is enough to consider only the zero modes and the first KK state. 



E.3.1 Incoming Spin-Down Massless Particle 

Thus we will assume that the incoming particle has energy slightly above (1/-R) and that it 
is reflected into an on shell first excited KK state while the on shell transmitted particle is 
massless. The 5D Dirac equation reads 



7° [V • 7 + V(x) + 7 5 d 5 ] 1> = Eij). 



(120) 



For a massless incoming particle with momentum k^ = (k, 0, 0, k) the solution is (again in the 
chiral basis) 
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Since the extra dimension is compact, there is no current flowing in this direction. The conserved 
current is therefore 



which implies, 



v • {r^) = o, 

R k \r\ 2 + \t\ 2 = 1. 



(122) 
(123) 



E.3.2 Incoming Spin-Down First KK Mode 



For a massive incoming particle with momentum k^ = (k, 0, 0, k') the solution this time is 
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As before, one obtains the relation 



\r\ 2 + R~ l \t\ 2 = 1. 



(124) 



(125) 



E.3.3 Incoming Spin-Up First KK Mode 

Here the solution is 



^{x 3 ) = 6{-x 3 ) 



Ak'x 3 



fVE=W\ 


y/WTk' 

V o / 



+ r\e 



—ikx 



l~3 



fVETw 
o 

y/E-k' 

V o / 



3\„ „-ikx 3 



+ 9(-x d )r e 



{V2E\ 



V o / 



J2 



and therefore the unitarity condition 



N 2 + tffeVol 2 = 1. 
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